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ABSTRACT 

There is a long-standing discrepancy between galaxy cluster masses determined from X-ray and 
gravitational lensing observations of which Abell 1689 is a well-studied example. In this work we 
take advantage of 180 ks of Chandra X-ray observations and a new weak gravitational study based 
on a Hubble Space Telescope mosaic covering the central 1.8 Mpc x 1.4 Mpc to eliminate the mass 
discrepancy. In contrast to earlier X-ray analyses where the very circular surface brightness has been 
inferred as Abell 1689 being spherically symmetric and in hydrostatic equilibrium, a hardness ratio 
map analysis reveals a regular and symmetric appearing main clump with a cool core plus some 
substructure in the North Eastern part of the cluster. The gravitational lensing mass model supports 
the interpretation of Abell 1689 being composed of a main clump, which is possibly a virialized cluster, 
plus some substructure. In order to avoid complications and mis-interpretations due to X-ray emission 
from the substructure, we exclude it from the mass reconstruction. Comparing X-ray and lensing mass 
profiles of the regular main part only, shows no significant discrepancy between the two methods and 
the obtained mass profiles are consistent over the full range where the mass can be reconstructed from 
X-rays (out to « IMpc). The obtained cluster mass within w 875 kpc derived from X-rays alone is 
6.4 ± 2.1 X 10^^ M0 compared to a weak lensing mass of 8.6 ± 3.0 x 10^** Mq within the same radius. 
Subject headings: galaxies: clusters: individual (Abell 1689), gravitational lensing. X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Gravitational lensing and X-ray observations are two 
independent methods, which can probe the mass of the 
largest virialized structures in the Universe, namely clus- 
ters of galaxies. The cluster mass is dominated by dark 
matter (« 80%) and hot X-ray emitting gas (~ 20%). 
The baryonic matter in the galaxies only contributes a 
few percent of the mass. The spectrum of the X-ray emis- 
sion from the gas depends on the temperature of the gas, 
which for a cluster in hydrostatic equilibrium is related to 
the total mass of the cluster. Gravitational lensing uses 
the distortion of background source images to probe the 
total mass along the line of sight. It has been widely 
debated whether cluster mass determinations from the 
two methods agree, a nd in some cases large discrepancies 



have b een foun d (e.g. Loeb fc Mao"1994| 'Miralda-Escude 
fc Babul 1 995 : Voi gt fc i''abTan 2006 : Mahdavi et al. 200 
Zhang et al. 2008). The main motivation for studying a 



i 



single cluster, Abell 1689, in large detail with the best 
available X-ray and lensing data is to get insight into the 
mass discrepancy by determining the mass distribution 
using both methods. It is checked that the results of a 
third independent mass determination method using the 
velocity dispersion of the galaxies within the cluster are 



danutas, 



desiree, 

marceau . limousiinaast .o bs-mip . f r , 



Electronic address: 

kp@dark-cosmology.dk, 

hakon.dahle@astro.uio.no< 

^ Dark Cosmology Centre, Niels Bohr Institute, University 
of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, 
Denmark 

^ Laboratoire d'Astrophysique de Toulouse- Tarbes, Universite 
de Toulouse, CNRS, 57 avenue d'Azereix, F-65000 Tarbes, France 

^ Institute of Theoretical Astrophysics, University of Oslo, P.O. 
Box 1029, Blindern, N-0315 Oslo, Norway 

Centre of Mathematics for Applications, University of Oslo, 
P.O. Box 1053 Blindern, N-0316 Oslo, Norway 



consistent with the X-ray and lensing results. 

Abell 16 89 is a massive c luster with a redshift of 
z = 0.1832 ("Teague et al.'1990'). It is mostly known for its 
amazing gravitational arcs and large number of multiply 
imaged systems. Consequently it has been well studied 
with strong and we ak gravitational lensing (jBroadhurstj 
et al.|2005a : Broad hurst et al |2005bl |Halkola et al. 200ff 



Lemze et al.||2007| |Limousin et al.|j2007( |Dahle et aT] 
2008). It has been proposed as a standard example of 
a spherical cluster of galaxies in hydr ostatic equilibrium 



dXue fc Wu||2002| [Lemze et al.||2007[ ). However, earlier 
studies have found large discrepancies between the mass 
obtaine d from X-ray obs ervations and f rom gravitational 
lensing (Andersson fc Madcjski 2004] [Andersson et ar 



2007 



Lemze et al. 2007). Strong gravitational lensing 



analyses have foun d the central 300 kpc to consist of 



several subclumps ( Limousin et al.||2007 Tu et al.||2007 l 

and weak lensing reveals substructure on a larger scale 
(500 kpc) (Okur a fc^ Futamasc 2008: Dahlc ct al. 2008j . 
Furthermore an Indication of substructure in the North 



Eastern (NE ) part of Abell 1689 was seen by Andersson 



fc Madejski (2004) in an XMM-Newton X-ray observa- 



tion. 

Recently new data have become available in the form 
of a very deep X-ray observation with Chandra (~ 150 ks 
new exposure). In addition, a weak gravitational lensing 
analysis based on a mosaic using the Wide Field Plan- 
etary Camera 2 (WFPC2) instrument aboard the Hub- 
ble Space Telescope (HST) has recently been carried out 



(Dahle et al. 2008). This provides a unique mass map 



with an unprecedented combination of spatial resolution 
(50") and a large spatial extension. We explore both of 
these data sets in this work. 

The long X-ray exposure allows us to create hardness 
ratio maps with a resolution of 5" (Sec. 2.4 1, revealing 
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TABLE 1 

The analysed Chandra observations of Abell 1689. 



Observation id 


Date 


Exposure time'' 


Data mode 


1663 


2001 Jan 7 


10.73 ks 


FAINT 


5004 


2004 Feb 28 


19.86 ks 


VFAINT 


6930 


2006 Mar 06 


75.79 ks 


VFAINT 


7289 


2006 Mar 09 


74.61 ks 


VFAINT 



^ The given exposure times are after light curve cleaning. 

substructure in the X-ray emission from the NE part of 
the cluster. In contrast, the South Western (SW) part 
seems circular and very regular. Consequently we as- 
sume it to be spherically symmetric and relaxed, and we 
reconstruct the temp erature and mass profiles from SW 
part alone (Sec. 2.5 1. The mass profile of the SW part 
is compar ed to the lensing mass profile of the same re- 
gion (Sec. |3.3| . The obtained profiles are in very good 
agreement and it is proposed that Abell 1689 consists 
of a spherically sym metri c part plus some substructure 

We conclude that the mass 
righ quality X-ray and gravita- 



in the NE part (Sec. 2.9 1 



profiles determined from 
tional lensing data are in agreement 
Throughout the paper we have used a standard ACDM 

cosmology with ^taryon = 0.046, = 0.72, ^CDM ^ 

0.23, and Ho = 70.1 kms~i Mpc'^ ( [Hinshaw et al 12008 1. 

At the redshift of Abell 1689, 1' corresponds to 185 kpc 
for the chosen cosmological model. Unless otherwise 
stated, quoted uncertainties are one a. 



2.1. 



2. X-RAYS 
X-ray observations 



We have analyzed four Chandra X-ray observations 
from the NASA HEASARC archive^ with a total expo- 
sure of approx. 180 ks (see Tab.[T|). It has been claimed 
that the two newest, and by far longest, observations 



canno t be used due to background issues (Lemze et al . 



20071, but these can be overcome, as shown in Sec. ^ 

The four observations in Tab. [l] have been reprocessed 
with the calibrati on CALDB 3 .4.1 and analysed with 
CIAO 3.4 (Fr uscione et al.|2006[ ) following standard pro- 
cedures^. We found the X-ray peak of Abell 1689 to have 
the position (R.A., decl.)=(197.87306°, -1.3413889°). In 
all observations, the X-ray peak was situated on the 
ACIS-I3 chip and we performed the analysis with events 
entirely from this chip. 



2.2. X-ray analysis: Image and surface brightness 

An exposure corrected X-ray image combined from the 
four observations is shown in Fig. [l] The emission is 
almost circular tempting one to conclude that Abell 1689 
is spherical and in hydrostatic equilibrium. However this 
is not the case, which is clearly demonstrated by the 
surface brightness profiles derived from Fig. [l] for the 
NE and SW halves respectively (defined by the dashed 
white line) . Fig. [2] shows the ratio between the NE and 
SW profiles. It is evident that within w 300 kpc the SW 
half is brighter, where outside « 500 kpc the two halves 
are more equally bright with a tendency for the NW half 
to be brighter. 
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Fig. 1. — Exposure corrected X-ray image composed from the 
observations in Tab. ^ (smoothed with a 10" wide Gaussian). The 
dashed white line through the X-ray peak (20° angle) divides the 
cluster in the NE and SW halves. The annuli described in Sec. 12.51 
are shown in solid green and the 3' circle in dashed green. 
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Fig. 2. — The ratio between the surface brightness profiles in the 
0.3 — 12.0 keV of the NE and SW halves with Poisson errors. 

2.3. X-ray analysis: Background 

Chandra has two telemetry modes called FAINT and 
VFAINT (Very FAINT). They have different event grad- 
i ngs and thereby different background even t selections 
(ICXO Proposers' Observatory Guide] |2007[ ch. 6.14). 
I'hree of the four observations analysed here have been 
done in the VFAINT mode. The standard blank sky ob- 
servations for background subtrac tion are produced from 
observations in the FAINT mode ( Markevitch|2003 1 , and 
hence the shape of the background spectrum from a 
VFAINT observation deviates from the spectrum of a 
blank sky observation. However, instead of using the 
pipeline processed VFAINT data, the data can be re- 
processed as FAINT observations. In the latter case, the 
shape of the background spectra of the two new, long ob- 
servations 6930 and 7289 are identical to the correspond- 
ing blank sky spectra. We have extracted the blank sky 
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spectra from the same chip regions as the observed clus- 
ter spectra and scaled the blank sky spectrum level to the 
corresponding observational spectrum in the 9-12 keV in- 
terval, where very little cluster emission is expected. 

Instead of using the blank sky background, it is pos- 
sible to use the local background from the same obser- 
vation allowing to take advantage of the improved back- 
ground reduction in the VFAINT mode. The background 
differs from chip to chip (in the CCD), so the local back- 
ground has to be extracted from the same chip as the 
source region (in this case ACIS-I3) . We have compared 
the two methods of background subtraction and the dif- 
ference between the final results (e.g. the temperature) 
is negligible. In the following we have used the blank sky 
method whenever extracting spectra, since it allows us 
to use more of the chip for the actual analysis and hence 
determine the temperature further from the centre of the 
cluster. For the images we have used the VFAINT data 
since the background event grading is better. 

2.4. X-ray analysis and results: Hardness ratio 

The ratio of low energy photons to high energy photons 
from the intra cluster medium is a proxy of the temper- 
ature structure of the cluster. If the cluster is isother- 
mal and only radiates thermally, the ratio between soft 
and hard photons (S/H) will be independent of position 
(within the cluster). However, if there is a variation in 
temperature, we expect to see a difference in S/H. Even 
clusters in hydrostatic equilibrium have a radial temper- 
ature gradient, which we expect to see as a circular struc- 
ture in the soft to hard photon ratio. 

From the two longest exposures (observation id 6930 
and 7289) we have created two hardness ratio maps with 
different energy splittings between the soft and hard pho- 
tons. One has a splitting energy of 1.0 keV so S/H = 
E[0.3-1.0 keV]/E[1.0-10.0 keV] and the other of 6.0 keV 
so S/H = E[0.3-6.0 keV]/E[6. 0-10.0 keV]. We produced a 
soft and a hard exposure corrected image with a binning 
of 5"/pixel. The images were smoothed with a Gaussian 
of width (T — 5", before the soft images were divided 
by the hard images. For visualisation, the hardness ra- 
tio maps were smoothed again with Gaussian of width 
a = 15", resulting in the two hardness ratio maps shown 
in Fig. [3] Bright colour means excess of low energy pho- 
tons. In the case of hydrostatic equilibrium this means 
colder, but since substructure is rarely in hydrostatic 
equilibrium, the temperature that can be derived from 
spectral fitting is not the actual physical temperature. 

In the 1.0 keV hardness ratio map (upper part of 
Fig. [3| a substructure in the NE part of the cluster is 
clearly visible, where in the 6.0 keV hardness ratio map 
(lower part of Fig. |3| only the (almost) circular symmet- 
ric part centered at the X-ray peak is visible. 



NE m 



m 



The split ting of 1 keV is commonly used (e.g. Fabian 
et al. 2000 1, but we have investigated several splittings 
between 0.5 keV and 7.0 keV. For splittings between 
1.0 keV and 6.0 keV, the hardness ratio map shifts grad- 
ually between the two extremes shown in Fig. [3] so in- 
creasing the splitting energy hides the NE structure and 
enhances the circular cluster structure. 

The contours of the total X-ray emission from 
Abell 1689 (green in Fig. [3]) give the impression of being 
circular with only a slight elongation in the North- South 
direction. This elongation is more pronounced in the 
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Fig. 3. — The obtained hardness ratio maps (bright is soft photon 
excess). The green contours show the almost circular total X-ray 
emission. The white line divides the cluster into the NE (upper) 
part and the SW (lower) regular and symmetrica lly a ppearing half. 
The white circle is the region analysed in Sec. |2.8| Upper: The 
hardness ratio map for the energy splitting of 1.0 keV (S/H = 
E[0. 3-1.0 keV]/E[1.0-10.0 keV]) Lower: The hardness ratio map 
for the energy splitting of 6.0 keV (S/H = E[0.3-6.0 keV]/E[1.0- 
10.0 keV]). 

hardness ratio maps with the 1 keV splitting. 

In the hardness ratio map with splitting at 6.0 keV 
there is a relative soft-photon excess centred at the posi- 
tion of the total X-ray peak, which indicates a cool core 
in the cluster. 

From the hardness ratio maps it is inferred that 
Abell 1689 consist of a spherically symmetric part, which 
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to the NE is obscured by come substructure. The cool 
core indicates that the spherical part of Abell 1689 is 
relaxed and in hydrostatic equihbrium (Voigt & Fabian 



2006 1 . Therefore we have split the cluster into two re- 



gions and in the following sections, the temperature and 
mass profiles have been determined from the symmetri- 
cally appearing SW part alone. The division line shown 
in Fig. |l] and Fig. [s] (white dashed) goes through the 
X-ray peak and is tilted by 20° with respect to the East- 
West direction. 

2.5. X-ray analysis and results: Temperature profile 

Assuming the intra cluster medium to be an opti- 
cally thin and completely ionised gas, we determined 
the global properties of Abell 1689 by analysing a 0.5 — 
8.0 keV spectrum of the central 3'. This radius was cho- 
sen because it includes most of the cluster emission (the 
dashed green circle in Fig. [I]) and the same radius is used 
in earlier an alyses ( Andersson fc Madejskij 2004| Lemze 
et al. '20071. An isothermal plasma model (MEKAL, 
Mewe et aT] (1985l) including Galactic absorption was 
fitted to the spe ctrum using the spe ctral fitting pack- 
age Xspec 12.3 (Arnaud et al. 20071. The absorption 



in neutral hydro gen along the l ine oi sight was fixed to 
1.83 X 102"cm-2 CHeaSoft«2007! uh tool). Leaving tlh a 
free parameter while fitting gave a consistent value. The 
redshift was fixed to z = 0.183. The redshift can also be 
determined from the spectral fitting, but the result had 
an uncertainty of « 20% and as a consequence we used 
the fixed value. Unfortunately this also excludes the idea 
of determining the distance along the line of sight to the 
NE substructure from the spectral fitting. 

We obtained a global temperature of 10.5±0.1 keV and 
an abundance of 0.37 ± 0.02 solar value for a reduced 
of 1.3 (1725 degrees of fr eedom). This is consistent with 
earlier reporte d results (Xue fc Wu 2002 Andersson & 
Madejs ki 2004- [Lemze etal.||2007| ). 

Clusters are not isothermal and consequently we made 
a radial temperature analysis. The X-ray emission of the 
NE part of Abell 1689 contain s sub structure in the low 
energies (see Sec. 2.4 and Sec. 
part of the cluster is complex. 



2.9 ) indicating that this 
Hence, we have restricted 
the radial temperature analysis to the SW symmetrically 
appearing half of the cluster. 

Using the four data sets and the blank sky back- 
ground, we have extracted the spectra from the SW half 
of Abell 1689 in 17 half-circular bins centered at the X- 
ray peak and containing at least 30000 events each (green 
circles in Fig.[T]). The same model as for the global prop- 
erties was fitted to the 0.5 — 8.0 keV spectrum of each 
radial bin. The obtained (projected, 2D) temperature 
profile is shown in Fig. |4] (black diamonds) . 

Fitting the complete annuli provides similar tempera- 
tures except for the annuli affected by the substructures. 
In general the reduced of the fits are closer to one for 
the half annuli than for the complete annuli. 

Since the observed X-ray emission is a superposition of 
gas shells of different temperatures we proceeded to de- 
project the temperature profile of Abell 1689 using ten 
half-circular bins, each containing at least 45000 counts. 
Using Xspec, the spectra and response matrices for all 
bins were fitted simultaneously to the same plasma model 
as for the 2D temperature profile, but now combined with 
the projct command, which accounts for the projection 



16 



14 



12 - 



fx 

i Tl' ■< 



2D projected 
3D deprojected 



8 - 



200 



400 600 800 
Radius [kpc] 



1000 1200 



Fig. 4. — The temperature profile reconstructed from the SW 
symmetric part of Abell 1689. The black diamonds are the 2D pro- 
jected profile and the grey crosses are the 3D deprojected profile. 
For the further analysis the 2D projected profiles was used. 

effects performing a 3D to 2D projection of prolate shell 
annuli. Each shell are assumed to have its own tem- 
perature and all the temperatures are allowed to vary 
independently. Spherical symmetry of the SW part was 
assumed as supported by the hardness ratio maps. The 
deprojected (3D) temperature profile is shown in Fig. |4] 
as grey crosses. 

The two temperature profiles are quite similar and for 
that reason, we have reconstructed the mass profile from 
the 2D projected temperature profile. 

2.6. X-ray analysis and results: Mass profile 

Assuming the SW part of Abell 1689 to be one half of 
a spherical symmetric cluster in hydrostatic equilibrium, 
the 2D projected temperatures can be used to determine 
a mass pr ofile following the procedure described in Voigt] 
& Fabian ( 2006) . The inethod is based on the hydrostatic 
mass equation ( Sarazin||l988 1: 



M3d{< r) 



hTg{r)r (d\n{Tg{r)) , d\n{pg{r)) 



Gfj^rUp 



c?ln(r) 



(iln(r) 



(1) 

where ks is the Boltzmann constant, G is the gravita- 
tional constant, /i is the mean particle weight weighted 
by the proton mass m^. Tg is the gas temperature as a 
function of radius, and pg is the gas density as function 
of radius, which can be found from the normalisation of 
the fitted plasma model. The resulting 3D mass profile 
of the SW part is shown in Fig. [5] (black) . The errors 
bars are determined by standard error propagation. 

2.7. X-ray analysis and results: NFW fit and projection 

The mass profile determined from X-rays describes the 
3D matter distribution. However the mass profiles from 
gravitational lensing (described in Sec. 3.11 are 2D pro- 
jected profiles, so in order to compare the two profiles, 
we need to project the X-ray profile as well. It is most 
easily done by projecting the density along line of sight, 
and t hen integrate along the projected radius, R (Sarazin 
T988|: 



M2d{<R)^ / 2^i?' 



Ptot{r)r 



drdU . (2) 
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Fig. 5. — The 3D X-ray mass profile (grey diamonds). The fitted 
NFW profile is shown as a grey dashed line, and the projected 2D 
NFW profile as a blue solid line. The 2D projected X-ray mass 
profile i shown as black diamonds. 

If we naively calculate the density profile from the dis- 
crete X-ray mass profile, it becomes unphysical with neg- 
ative densities because of the fluctuations in the mass 
(even if the uncertainties are taken into account). In- 
stead we have investigated how large the projection effect 
is for different analytical density profiles. The projection 
mainly affects the inner part of the mass profile by a 
factor of up to two, and leave the outer parts effectively 
unchanged. In the very central parts of the cluster, the 
projection depends on the cuspiness of the chosen profile. 
However, the X-ray mass profile does not probe the cen- 
tral part in great detail and the chosen profile is required 
to reproduce the observed 3D mass profil e by fitting, so 



we h ave chosen to use an NFW profile (Navarro et al 
19961. The 3D NFW integrated mass profile was fitted 
to the 3D X-ray mass profile using reduced statis- 
tics. The best fit is shown as a dashed grey line in Fig. |5] 
and has the following parameters: = 174 ± 10 kpc, 
po = (7.79 ± 0.02) X lO^Mo/kpc^, csoo = 5.6 for a re- 
duced = 1-6 (13 de gr ees of freedom) con sistent with 
earlier pubhshed results ( [Lemze e t "ar] l2007[ ). The NFW 
profile was projected using Eqn. [5] and the obtained 2D 
profile was divided by the 3D profile to achieve a radius- 
dependent projection factor, which ranged from approx. 
two in the centre of the cluster to one in the the outer 
parts. Each point in the 3D X-ray mass profile was con- 
verted to a 2D mass profile by multiplying with the pro- 
jection factor at the corresponding radius. Fig. |5] shows 
the 3D profile as grey diamonds, the fitted NFW pro- 
file as a grey dashed line, the projected 2D NFW profile 
as a blue solid line, and the corrected 2D X-ray mass 
profile as black diamonds. The one a error bars were 
correspondingly corrected. 

2.8. X-ray analysis and results: The NE region 
The NE region is dominated by soft photons (see 



Sec. 2.4 1 and for that reason expected to be colder than 
the SW region. Fitting a single temperature component 
plasma model (same model as above) to the spectrum of 
the region within the white circle in Fig.[3]gives a temper- 
ature of 9.3 ± 0.9 keV (reduced = 0.97 for 129 d.o.f.). 
This is slightly colder than the corresponding SW part 
of the cluster, which has a temperature of 11.1 ± 0.6 keV 
(reduced = 1-02 for 707 d.o.f.) at similar distance to 



the centre. The emission from the region is assumed to 
be a combination of the emission from the SW spherical 
component plus some softer component. 

We have attempted to determine the temperature of 
the NE component, by fitting a two temperature model 
to the spectrum of the NE region. One of the tempera- 
tures was fixed to the temperature of the spherical com- 
ponent at the same distance (11.1 keV) while the other 
temperature was a free parameter of the fit. We ob- 
tained a cold component temperature of 1.3 ± 1.0 keV 
(reduced = 0-92 for 128 d.o.f.) with a normalisation 
of 9% of the normalisation of the warm component. Maz- 



zotta et al. ( 2004 1 demonstrates that in general it is very 
hard to disentangle two temperature components with 
equal normalisation unless in the case where one tem- 
perature is around 10 keV and the other one well below 
a few keV. Comparing to our scenario, the normalisations 
are not equal. Therefore we have tested the stability of 
the two temperature fit by varying the two temperatures 
and their fractional normalisation and then fitting for 
the overall normalisation. This procedure gives similar 
reduced for a range of parameters. For a reduced 
X^ < 1-5 the two temperatures can be varied indepen- 
dently between 4 keV and 14 keV as long as the normal- 
isations stays a free parameter. The conclusion is that 
the best fit for the NE region is either a two tempera- 
ture model with a warm and a very cold component or a 
one temperature model with a temperature slightly lower 
than the corresponding SW part of the cluster, which is 
consistent with fittin g a one temperature m odel to a two 
temperature plasma ( Mazzotta et al.|[2004 l. 



2.9. Discussion of X-ray results 

The hardness ratio maps and spectral analysis reveal 
that Abell 1689 is not spherically symmetric as indicated 
by the total X-ray image, but rather consists of a spher- 
ical main clump centred at the X-ray peak and a softer 
X-ray emitting substructure to the NE. 

The NE substructure is not easily visible in the total 
X-ray emission. It is because the total cluster X-ray emis- 
sion is generally dominated by low energy photons almost 
regardless of gas temperature, which again is dominated 
by emission from the main cluster. The absolute number 
of photons leading to the relative difference in the soft to 
hard photon ratio is simply drowned in the total X-ray 
emission. However, the substructure is clearly visible in 
the ratio between the surface brightness profiles of the 
two halves shown in Fig. [2] 

The substructure visible in the hardness ratio in the 
upper part of Abell 1689 is consist ent with the results 
from XMM-Newton presented by Ander sson fc Madejski] 
p004), but it is more distinct in the work presented here. 
T'he symmetric SW part of Abell 1689 has a cool core 
indicated by the hardness ratio map and by the tem- 
perature profile. The cool core has earlier been cl aimed 

2004| . The 
TTB89 is in 



missing in Abell 1689 (Andersson & Madejski 
cool core indicates, that the SW part of Abe 



hydrostatic equilibrium ( Voigt &: Fabian||2006 1 . 

The temperature profile obtained from the SW part of 
the cluster is compared to earlier published X-ray tem- 
perature profiles in Fig. [61 The error bars are smaller 
than earlier obtaine d profiles from Chandra data ( Xue & 
Wu|2002 Lemze et al.,2007) due to better statistics, and 
the measurements extend to larger radii due to the use 
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Fig. 6. — The obtained tem perature profile co mpared to earlier 
results: purple (dashed) from'Xue fc Wu'( j2002[ l, orange (dashed 
dotted) from Andersson & Madejski (2UU4) (ijL) deprojected), and 
red (dotted) from , Lemze et ai.( ( |2007J . 

of blank sky ba ckground. Only the an alysis of XMM- 
Newton data by Andersson & Madejski (2004) has com- 
parable statistics and extension. However the temper- 
ature profile obtained in this work is generally slightly 
warmer which leads to an overall larger mass. 

The X-ray mass profile as inferred from the SW part 
of the cluster is well fitted by an NFW profile with a 
concentration parameter of C200 = 5.6. This values is 
in better agr eement with the value of C200 = 4.7 ±1.2 

(20071 from cosmological simu- 
7.6 ±1.6 obtained from 
I or the value 

of C2nn = 7.7 fr om earlier X-ray analyses Andersson Sz\ 



predicted by Neto et al 
lations, than the value of C200 



gravitational lensing ( Limousin et ar|[2007 
7.7 from earlier X-ray analyses 



Madejski (20041. This supports our interpretation of the 
SW part bemg spherical and in hydrostatic equilibrium. 

3. ADDITIONAL OBSERVATIONS OF ABELL 1689 
3.1. Strong and weak lensing analyses 

We compare our X-ray analysis to strong and weak 
gravitational lensing analyses of Abell 1689. Strong and 
weak gravitational lensing are two different methods of 
determining the mass of the lensing object (s) based on 
distortions of the images of background sources. Strong 
lensing is tracing the mass of the lens within the Ein- 
stein radius where arcs and multiple images arc forming. 
Weak lensing traces the mass outside the Einstein radius 
where the gravitational force is weaker, and stretches the 
images of background objects perpendicular to the di- 
rection to the center of the lens. Thus the two methods 
are valid in different ranges, overlapping only around the 
Einstein radius. 

There ha ve been many strong lensing studies of 



Abell 1689 QBroadhurst et al.||2005a| [Broadhurst et al. 
2005b| [Limousin et al.|2007| [Lemze et al.|2007p. Here we 



have u sed the latest and most detailed by |Limousin et al.[ 
(20071, which is based on data from the HST Advanced 
Camera for Surveys (ACS) with the spectroscopic infor- 
mation of the lensed objects in the system from the Keck 
Telescope and the Very Large Telescope (VLT). The 
mass reconstruction was done using a Bayesian Monte 
Carlo Mark ov Chain method with LENSTOOL (Tdlo 
et al.[[2007l. The result is an accurate mass mapping 




Fig. 7. — Optical images (HST/ACS) with gravitational lensing 
contours. The black contours are the mass contours from strong 
lensing. The minimum value is 10^ Mq/"^ and the distance be- 
tween contours is 3.5 X 10* Mq/"^. The strong lensing has a mass 
peak at the centre of the cluster and a smaller peak ^ 180 kpc to 
the NE (A). The white contours are the k contours from weak grav- 
itational lensing with a resolution of 50" (FWHM). The minimum 
value is K = 0.1 and the distance between contours is Sk = 0.16. 
The weak lensing contours clearly shows a subclump 500 kpc to 
the NE of centre (B). 

a clear bimodality as shown in Fig. [T] (black contours). 
In the very centre there is a main clump associated with 
the peak of the X-ray emission and then a second signif- 
icant clump « 180 kpc NE of the centre (A in Fig. [tI 
associated with a clump of galaxies (not to be confused 
with the X-ray structure to the NE, at a di stance from 
~ 500 kpc to the centre). Previo us studies (Broadhurst 
eral]|2005b[ [Halkola et aL]p006] ) were using fewer spec- 
troscopically confirmed imaged background sy stems but 
are in good agreement with the results from Limousin 



of the central parts of Abell 1689, which demonstrates 



et al.[ ( [2007 l. 

Tot weak grav itational lens i ng we compare to a re- 
cent analysis by [Dahle et al.[ (20081, which was based 
on a mosaic of 16 HST WFT'C2 pointings covering the 
central ~ 1.8 Mpc x 1.4 Mpc of the cluster, comple- 
mented at larger radii by ground-based data using the 
CFH12K mosaic CCD camera at the 3.6 m Canada- 
France-Hawaii Telescope (CFHT), covering a larger field 
of 7.7 Mpc X 5.2 Mpc. This provides a unique mass map 
with an unprecedented combination of spatial resolution 
(50") and a large spatial extension. The surface mass 
density map from weak gravitationa l lensing shown in 
Fig. [t] (w hite) was reconstructed by [Dahle et al. ( 2008 1 
using the|K aiser & Squires[ ([1993[) method. This method 
relies on the fact that the convergence and the two com- 
ponents of the shear are linear combinations of the sec- 
ond derivative of the effective lensing potential. Using 
Fourier transformation of the convergence and the shear, 
one can obtain linear relations between the transformed 
components and determine the mass. The mass contours 
clearly shows a non-spherical morphology of Abell 1689 
at « 500 kpc from the centre (B in Fig. [8|. 
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Fig. 8. — The hardness ratio maps with energy spHtting of 
1.0 keV (S/H = E[0.3-1.0 keV]/E[l. 0-10.0 keV]). Bright is soft pho- 
ton excess. The white Hne divides the cluster into the NE (upper) 
part and the SW (lower) symmetric part. The contours are identi- 
cal to Fig. [7] with white being weak gravitational lensing and black 
strong gravitational lensing. 

3.2. Lensing analysis: mass profiles 

From weak lensing the one dimensional radial mass 
rofile was obtained using aperture mass densitometr y 



Kaiser 1995 Clowe et al. 2000 Fahlman et al. 19941 



his method determines the mean surface mass density 
within an aperture minus the mean density in a sur- 
rounding annulus. Hence, a lower limit to the projected 
mass of the lens inside the aperture can be found. The 
degree of underestimation of the true enclosed mass de- 
pends on the inner and outer radii and the slope of the 
projected mass density profile. For realistic projected 
density profiles (steeper than at large radii) the mass 
results presented here will be underestimated by less than 
30% at all radii. 

For strong gravitational lensing the mass map was in- 
tegrated in radial bins of 1" centered on the peak of the 
X-ray emission (which coincides with the strong lensing 
mass peak). 

For both weak and strong lensing we reconstructed 
two different radial mass profiles: one by integrating the 
whole surface density map and the other by integrating 
the SW half of the surface density map as in the X-ray 
analysis (but multiplied by two to get a complete sphere) . 

3.3. Lensing results and discussion 

Both strong and weak lensing analyses show that 
Abell 1689 is not spherical, but has substructure in the 
NE direction on both small (w 180 kpc) and large scales 
(« 500 kpc). 

The mass profiles of the SW part of Abell 1689 ob- 
tained from gravitational lensing are shown in Fig. 9] 
The mass profiles were plotted keeping in mind that the 
strong lensing mass map is valid within the Einstein ra- 



Lensing: Total cluster 
Lensing: 2 times SW half 
X-ray main clump 2D 




100 



1000 
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Fig. 9. — The 2D projected mass profile reconstructed from X-ray 
(black), strong gravitational lensing (blue) and weak gravitational 
lensing (green). The dashed lensing profile is for integration over 
the full cluster and the lower solid curve is for integration over the 
SW part multiplied by two (in order to get a full sphere). The 
width of the shaded regions corresponds to the one sigma errors on 
the lensing profiles. 

dius, which for Abell 1689 is 45 139 kpc. The dashed 
curves in Fig. |9] are the integration of the whole mass 
map and the solid curves are the integration of the SW 
part only, but multiplied by a factor of two to get a full 
sphere. 

For the strong lensing, the mass profile integrated over 
the total cluster is very similar to the mass obtained by 
integrating the SW part only since the NE substructure 
at « 180 kpc is negligible compared to the total mass 
of the cluster. For the weak lensing, the difference is a 
bit smaller and the total integration gives a larger value 
due to the NE substructure also seen in the X-rays. The 
total mass within 1200 kpc is « 9 x lO^"* Mq. 

The mass of the weak lensing substructure is deter- 
mined to be 1.25 ± 0.3 x IO^-^Mq within 1.56 '. 

In Fig. [To] the total cluster mass profile is compared 
to earlier studies of Abell 1689 with weak gravitational 
lensing. It is seen that there is good agreement between 
the strong lensing mass profiles for radii smaller than 
100 kpc. Above 100 kpc, the profile of Dahlc et al. 
(2008) is slightly lower than the earher profiles ([Broad 



hurst et al 2005a, ,Halkola et al.|2007| [Umetsu fe Broad- 



hurst:2007l 



Assuming the strong lensing mass profile to be the 
actual 2D projected mass profile of the inner part of 
Abell 1689 and the X-ray profile to be the actual 3D pro- 
file allows us to constrain the 3D to 2D projection factor 
as a function of radius, which depends on the cuspiness 
of the mass profile. The maximally allowed projection 
factors are M2D/M3D = (2.8,4.2,3.1,2.6) for radii of 
r = (38, 67, 92, 188) kpc, which are all larger than the 
factors determined from the NFW profile projection in 
SeclTzl 

3.4. Redshift distribution of cluster galaxies 
The redshift distribution of the galaxies within 



Abell 1689 is ov erall Gaussian (Due et al.|2002 Teague 



et al.||1990 [[Cibskc 2002), but witF irregularities (|Lokas| 
et al. S()06). The full width at half maximum of the mam 
structure is Az = 0.02 corresponding to a line of sight 
velocity dispersion of w 2500 km/s. This is in contrast 
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Fig. 10. — The weak and strong gravitational mass profiles 
(dashed navy) of t he tota l cluster compared to earlier results by 



[Broadhurst et al.| | |2005a^ (oliye green) and Halkola ct al. (2007) 
(green) . Also the A-ray rnass profile obtain ed here is shown ( black 
crosses) together with the one obtained by [Andersson fc Madejski] 
|2004| (aquamarine dashed). 

to other clusters where very few rich clusters have been 
found with a velocity dispersion above 1000 — 1200 km/s. 
It suggests either that Abell 1689 is not in hydrostatic 
equilibrium or it is a structure with several clumps along 
th e line of sight. The latter i nterpretation is favo ured 



by Girardi et al 
compared the ve 



( |1997[ ) and by |Lokas et al.] |2006[ | who 
ocity distribution of Abell 1689 to sim- 
ulations and argued that the high velocity dispersion is 
due to several clumps along the line of sight. Remov- 
ing what they found to be separate unbound structures, 
they got a total dynamical mass of 2 x 10^^ M© within 
3.5 Mpc. 



Also Lemze et al. ( 2008 ) have studied the dynamics of 
Abell 1689. They derived the virial mass of the cluster 
to be 1.73 ±0.59 x 10^^ M© using the velocity caustics of 
the cluster galaxies. 

4. DISCUSSION 

Comparing the X-ray hardness ratio maps with the 
weak and strong gravitational lensing mass maps (Fig.^ 
we see a coincidence between the substructure in the X- 
ray emission and the bimodality in the lensing maps. The 
total X-ray emission (Fig. [T|) features a slight elongation 
along the same axis as the secondary clump in the strong 
lensing map. Further from the centre, the position of the 
NE substructure in the hardness ratio map coincides with 
the position of a large substructure in the weak lensing 
map. On the opposite site, the SW part of the cluster 
appears circular and symmetric in all observations. From 
this we conclude that Abell 1689 consist of a spherical 
main clump and some substructure to the NE, which is 
interfering with the emission from the NE half of the 
main clump. 

The NE half of the cluster is slightly brighter than the 
SW half at the distance of the weak lensing substructure 

The mass of the substructure is 1.25 ± 0.3 x IOI-^Mq 
within a radius of 1.56 'determined from weak lensing. 
This is only a small fraction of the total cluster mass. 
The gas mass within the same radius has been estimated 
assuming sphericity and uniform gas density. The result 
of sa 9.7 ± 0.3 X 10^^ M© only contributes with a very 



small fraction to the total cluster mass. Nonetheless the 
X-ray emission from the substructure is bright enough 
to introduce a complication of the spectral analysis and 
thereby potentially introduce a mis-interpretation of the 
temperature of the cluster. 

Reconstructing mass profiles from the SW half of the 
main clump only leads to unprecedentedly good agree- 
ment between weak and strong gravitational lensing, and 
the 2D projected X-ray mass profile, which supports 
the interpretation of Abell 1689 consisting of a spheri- 
cal main clump and substructure to the NE. 

We have quantified the deviation between the differ- 
ent mass profiles by comparing the relevant part of the 
2D projected X-ray mass profile to the weak and strong 
lensing profiles separately. The values of one profile were 
compared to a linear interpolation of the other profile. 
The reduced was determined to be 1.8 (four degrees 
of freedom) for the weak lensing and X-ray comparison 
and 1.2 (10 degrees of freedom) for comparison of the 
strong lensing to the X-ray mass profile. 

The cool core indicates that Abell 1689 is not a re- 
cent merger and supports the claim that the SW main 
clump is in hydrostatic equilibrium. It indicates that the 
NE substructure is not leftovers from an early collision 
and maybe not gravitationally interacting with the main 
clump at all. If it is infalling, the merging is in such 
an early state that the hydrostatic equilibrium of the 
SW part has not yet been disturbed. The line of sight 
galaxy velocity distribution supports a scenario with sev- 
eral clumps along the line of sight, where the gas does 
not interact. 

The hardness ratio maps can be used as easy diag- 
nostics for the thermal distribution of the cluster, which 
is related to any substructure. For Abell 1689 the sub- 
structure to the NE is visible in a hardness ratio map 
with an energy splitting of 1.0 keV and a resolution of 
10" made from 15 ks of the 6930 observation. Many clus- 
ter observations have exposures of the order of 15 ks so 
we propose to use hardness ratio maps in the selection 
process of relaxed clusters to be further studied. 

In this work we have taken advantage of high quality X- 
ray and gravitational lensing data, which have improved 
the agreement between several mass estimation methods 
significantly. 

5. CONCLUSIONS 

The new X-ray data analysed here and compared with 
weak and strong gravitational lensing data show that on 
large scale Abell 1689 appears spherical, in agreement 
with earlier data. However, both temperature structure 
and gravitational lensing shows that Abell 1689 contains 
some substructure to the NE. The main part appears cir- 
cular and is centered at the peak of the total X-ray emis- 
sion. Only the SW half is not significantly influenced by 
other structures. From X-ray hardness ratio maps and 
temperature profile it is seen that the SW part features a 
cool core. The NE substructure is seen in the weak grav- 
itational lensing mass map and the X-ray hardness ratio 
maps. It has an excess emission of soft X-ray photons 
relative to hard photons. 

We have determined the mass profile of the spherical 
main clump of the cluster from X-ray observations and 
compared to recent gravitational lensing results. The ob- 
tained cluster mass within w 875 kpc derived from X-rays 
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alone is 6.4 ± 2.1 x 10^^ Mq compared to a weak lensing 
mass of 8.6 ± 3.0 X 10" Mq within the same radius. The 
profiles are in very good agreement out to 1200 kpc and 
the discrepancy between X-ray and lensing mass profiles 
has been significantly reduced due to high quality data. 
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